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We show that the dominant channel proposed for the determination of the Higgs boson trilinear 
coupling, pp —> HH + X via gluon fusion, exhibits an interference structure that is independent of 
the collider energy for collider energies in the range 8 TeV < y/s < 100 TeV and is almost maximally 
destructive. This insensitivity to the collider energy remains approximately true for a variety of 
other two Higgs production mechanisms although the magnitude of the interference varies widely. 

PACS numbers: 13.38.Dg 


1. INTRODUCTION 

The standard Higgs self interaction is contained in the Higgs potential 

V{H) = Xv 2 H 2 + XvH 3 + ^XH 4 , (1) 

where A = m 2 H /2v 2 and, in terms of the W mass, the weak mixing angle and the fine structure constant, 
v = Mw sin Ow/y/na. A measurement of the trilinear coupling is an important test of the Standard model 
behavior of the Higgs-like object discovered at the LHC [1, 2]. Given the proliferation of gluon production at the 
Large Hadron Collider, the natural choice for the study of the Higgs boson trilinear coupling is the gluon fusion 
process gg -A HH [3-13]. The matrix element comes from the sum of a triangle graph and several box graphs 
as indicated in Fig. 1. Each of these two contributions is gauge invariant. By calculating the total cross section, 



FIG. 1: The leading order triangle and a representative box diagram are shown. 

c ttot and separately calculating the cross sections from just triangle graph, <tt, and just the box graphs, os, 
we can study the interference between the two amplitudes. Defining interference angle ai by the equation 

otot — + 0's + 2 cos (aj) y/arap , (2) 

we can examine how cos(aj) varies as a function of energy. Since the interference term between the box and 
triangle amplitudes is 2 bM = 2 \M. B \\M. t\ cos(a/), cos («/) is independent of the three Higgs boson 
self coupling. 
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2. RESULTS FOR gg -> HH 


The numerical results that are shown in Table I were obtained using Mh = 125.5 GeV, M t = 173.1 GeV and 
the MSTW2008 [14] parton distribution functions (PDFs) with both scales equal to the invariant mass of the two 
Higgs final state. The leading-order amplitudes for Higgs pair production from gluon fusion can be expressed in 
terms of tensor integrals [15] and scalar integrals [16]. We evaluate these integrals numerically with a FORTRAN 
code [17] developed for this purpose. The cross section was calculated in the usual way and then Eq. (2) was 
used to find cos(aj). The error in these values of cos (a/) from the numerical integration varies from 0.005 at 
low energies to 0.01 at high energies. The negative value of cos(o:j) reflects the known fact that the interference 


Vi(TeV) 

rj B (fb) 

a T (fb) 

CTTOT(fb) 

cos(ar) 

8 

9.06 

1.34 

4.11 

-0.902 

13 

31.6 

4.36 

14.9 

-0.898 

14 

37.8 

5.16 

17.9 

-0.898 

33 

243. 

30.3 

120. 

-0.893 

60 

760. 

89.8 

383. 

-0.893 

100 

1900. 

212. 

965. 

-0.904 


TABLE I: For pp —> HH + X from gg —> HH, the center of mass energy, the LO contributions to the cross section from 
the box, the triangle, and the total cross section, and cos(au). 

is destructive. What is surprising is that cos(aj) is almost universal. 

To check the numbers in Table I, we multiply the total cross section, otot, by the K factor given by de 
Florian and Mazzitelli [12] 

/ i? \ — 1 / e \ _1 / 2 

A ' = L242 - 7 - 17 (r5v) + 5,77 (livv) ■ (3) 

where E is the center of mass energy. This reproduces their NNLO numbers precisely. 

Since it is not clear how cos (ay) can be so constant with energy, it makes sense to ask how it varies with 
other factors. It is well known that the process, gg —> HH is very sensitive to the values of the scales and the 
PDFs. To see how cos(a/) depends on the scales and PDF we give in Table II the box contribution, triangle 
contribution, total cross section, and cos (aj) for three combinations of scale and PDF at three energies. For 
each energy the first line is the MSTW PDF with the scales equal to the two Higgs invariant mass (same as 
Table I). The second line uses the CTEQ6L1 PDF [18] with the same scales, and the third line uses the MSTW 
PDF but with the scales equal to Mh- Again the cross sections were calculated from the matrix elements and 
then cos(a/) was determined from Eq. (2). 


y/s (TeV) 

ffs(/6) 

a T (fb) 

ctot (/ 5 ) 

cos (a/) 

8 

9.06 

1.34 

4.11 

-0.903 


8.05 

1.20 

3.64 

-0.902 


15.3 

2.06 

7.32 

-0.894 

14 

37.8 

5.16 

17.9 

-0.897 


35.0 

4.79 

16.5 

-0.899 


60.5 

7.53 

29.9 

-0.893 

100 

1900. 

212. 

965. 

-0.904 


1790. 

198. 

914. 

-0.902 


2530. 

266. 

1340. 

-0.887 


TABLE II: For the cross sections of pp —> HH + X from gg —» HH the dependence of cos (ai) on the PDF (first and 
second line at each energy) or on the scale (first and third line at each energy). 

These results are very striking because, as can be seen, the contributions to the triangle and box cross sections 
and the total cross section are very dependent on the scales and somewhat dependent on the PDF. But cos (or) 
apparently doesn’t care about such things. 

To separate the effects of the underlying physics from that of the parton distribution functions, consider the 
contributions with the distribution functions set to unity, effectively examining the behavior of the partonic 
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process. Here we have to consider smaller energies because the triangle graph gets small rapidly due to the 
s-channel pole. The results are shown in Table III. 


VS (GeV) 

VB(fb) 

a T (fb) 

OTOt(/&) 

cos(au) 

300 

0.199 

0.682X10" 1 

0.342X10- 1 

-1.00 

400 

0.964 

0.875x 10 _1 

0.482 

-0.980 

500 

0.900 

0.497X10" 1 

0.552 

-0.940 

1000 

0.199 

0.277x 10~ 2 

0.165 

-0.783 

2000 

0.447X10” 1 

0.995x 10~ 4 

0.425X10" 1 

-0.545 

4000 

O.lMxlO" 1 

0.303x 10 -5 

0.113x 10 _1 

-0.277 


TABLE III: The cross sections and cos (a/) with the distribution functions are set to one, i.e., the parton cross sections 
a(gg HH) are shown. 

The last two lines are uncertain because they depend on the <jb and utot being very accurate. However, the 
conclusion is clear; without the distribution functions to emphasize the low energy parts of the cross sections, 
cos(oy) is no longer universal. More specifically, this Table shows that the behavior of the partonic cross section 
near threshold (VS ~ 300 — 500 GeV) deduced in [19, 20] using unitarity cutting rules is the basically the only 
region probed by the full cross section a (pp —> HH + X) no matter what the center of mass energy happens to 
be. 


3. INTERFERENCE IN OTHER TWO HIGGS PROCESSES 

3.1. Production of two Higgs in association with a Z or W 


It is worth looking at other two Higgs production processes to see how cos(a/) varies with energy for them. 
We call the background, B, the contribution from the diagrams where both Higgs are emitted separately from 
the Z line. The signal is the contribution from the diagram with the three Higgs coupling. In analogy to the 
gluon fusion case above, we call that T. The total cross section is the square of all the diagrams including the 
cross terms between the B part and the T part, cos(ay) is then defined by Eq. (2) above. First consider the 
results for pp —> ZHH + X shown in Table IV for the same energies used in Table I. 


\fs (TeV) 

&B 

(T T 

utot 

cos (or) 

8 

5.37x 10~ 2 

9.10x 10~ 3 

9.38x 10~ 2 

0.701 

13 

0.135 

2.21x 10~ 2 

0.229 

0.658 

14 

0.153 

2.51xl0” 2 

0.260 

0.661 

33 

0.588 

9.31x 10~ 2 

0.975 

0.628 

60 

1.34 

0.212 

2.23 

0.636 

100 

2.67 

0.411 

4.34 

0.601 


TABLE IV: Cross sections and interference for pp —¥ ZHH + X. The cross sections are in fb. The different contributions 
again refer to Eq. (2) above. The error in cos(aj) due to roundoff is less than or equal to 0.01. 

We see that cos (a/) changes only from 0.7 to 0.6 so unlike gluon production of two Higgs it does change but it 
doesn’t change much. For this reaction and the remaining reactions, the CTEQ6L1 distribution functions were 
used with scale equal to V§- We did try varying the distribution functions and found that it makes very little 
difference which set of PDF’s is used. The scale used makes large differences in the values of cross sections but, 
for a given energy, they all change in the same way so cos (aj) is unchanged. 

The processes pp —> W^HH + X have an interference behavior similar to pp — > ZHH + X. For the same 
range of energies, cos(oy) changes from 0.66 at 8 TeV to 0.57 at 100 TeV for W + and 0.69 to 0.58 for W~. The 
results are shown in Tables V and VI. 

For pp — > ZHH , rather than pp — > ZHH + X as given above, cos (a/), shown in Table VII, is 0.63 at 8 TeV 
and 0.60 at 100 TeV. So the change in the interference with energy when using two quark distributions is even 
less than when using one quark and one antiquark. 














4 


y/s (TeV) 

&B 

(Tt 

(JTOT 

cos ( q /) 

8 

5.81x 10~ 2 

1.15xl0 -2 

0.104 

0.665 

13 

0.136 

2.61x 10~ 2 

0.237 

0.629 

14 

0.153 

2.92x 10” 2 

0.266 

0.627 

33 

0.543 

0.101 

0.923 

0.596 

60 

1.21 

0.221 

2.03 

0.579 

100 

2.31 

0.418 

3.85 

0.571 


TABLE V: Cross sections and interference for pp —» W + HH + X. The cross sections are in fb. 


yfs (TeV) 

& B 

CJ T 

otot 

cos (or) 

8 

2.54x 10~ 2 

5.10x 10“ 3 

4.63x 10~ 2 

0.694 

13 

6.84x 10~ 2 

1.35x 10~ 2 

0.122 

0.660 

14 

7.86xl0~ 2 

1.54xl0” 2 

0.139 

0.647 

33 

0.333 

6.28x 10~ 2 

0.573 

0.613 

60 

0.807 

0.150 

1.37 

0.594 

100 

1.64 

0.300 

2.76 

0.585 


TABLE VI: Cross sections and interference for pp —> W HH + X. The cross sections are in fb. 


yfs (TeV) 

°B 

ctt 

otot 

cos(au) 

8 

0.104 

1.66xl0~ 2 

0.173 

0.631 

14 

0.217 

3.40xl0~ 2 

0.360 

0.634 

33 

0.657 

0.102 

1.07 

0.601 

60 

1.42 

0.219 

2.31 

0.602 

100 

2.72 

0.417 

4.41 

0.598 


TABLE VII: Cross section and interference for pp —> ZHH. The cross sections are in fb. 


3.2. Production of two Higgs in association with a tt pair 


Another reaction for producing two Higgs is pp —> ttHH + X. This has contributions from initial quarks and 
from initial gluons which can not be separated experimentally. But they must be calculated separately so let 
us first consider the qq —> ttHH part. For all the same input parameters, distribution functions, and scales as 
above our results are shown in Table VIII. Here, cos(a/) varies only from 0.84 to 0.76 over the usual enormous 
energy range. The other contribution to pp -A ttHH + X is from gg —> ttHH. The results are shown in Table 
IX. Of course we can’t measure the quark and gluon processes separately. If we ignore the fact that the K factors 


Ti(TeV) 

&B 

ctt 

otot 

cos (or) 

8 

6.47x 10~ 2 

1.96x 10~ 3 

8.56X10' 2 

0.841 

13 

0.212 

6.54xl0” 3 

0.280 

0.825 

14 

0.249 

7.68 x 10~ 3 

0.328 

0.815 

33 

1.19 

3.75X10" 2 

1.56 

0.788 

60 

2.98 

9.51xl0~ 2 

3.89 

0.765 

100 

6.15 

0.197 

8.04 

0.769 


TABLE VIII: Contributions to the cross section pp —> ttHH + X from qq —> ttHH. As always the cross sections are in 
fb. 

would be different and just add these last two LO processes to get a total cross section for pp —► ttHH + X we 
get the result shown in Table X. Finally we have a process where cos(ay) varies substantially with the energy. 
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y/s (TeV) 

&B 

<7T 

otot 

cos(aj) 

8 

2.79x 10 -2 

1.40xl0” 3 

3.28x 10“ 2 

0.280 

13 

0.179 

1.03x 10~ 2 

0.211 

0.253 

14 

0.231 

1.36xl0” 2 

0.273 

0.253 

33 

2.93 

0.205 

3.47 

0.216 

60 

12.6 

0.966 

15.0 

0.206 

100 

38.2 

3.11 

45.5 

0.192 


TABLE IX: Contributions to the cross section pp —> ttHH + X from initial gluons. 


Ti(TeV) 

c tb 

<JT 

otot 

cos (a/) 

8 

9.26x 10 -2 

3.36x 10 -3 

0.120 

0.681 

13 

0.391 

1.68xl0 -2 

0.491 

0.513 

14 

0.480 

2.12xl0' 2 

0.600 

0.490 

33 

4.12 

0.243 

5.03 

0.333 

60 

15.6 

1.06 

18.9 

0.275 

100 

44.4 

3.31 

53.6 

0.243 


TABLE X: Contributions to the cross sections and interference for the total process pp —> ttHH + X. 


3.3. HII production from qq —> qqHH 


The results for pp —> qqHH via the subprocesses uu uuHH, dd —> ddHH and ud udHH are shown 
in Tables XI,XII and XIII. Again, there is very little spread the values of cos (a/) but the interference is, like 
gg -A HH, destructive. 


Ti(TeV) 

OB 

CTt 

(ttot 

cos (cm) 

8 

5.87x 10~ 2 

1.53x 10~ 2 

2.30x 10~ 2 

-0.846 

13 

0.161 

3.81x 10“ 2 

6.90x 10~ 2 

-0.830 

14 

0.185 

4.32x 10 -2 

8.43xl0' 2 

-0.805 

33 

0.766 

0.157 

0.376 

-0.789 

60 

1.80 

0.341 

0.938 

-0.768 

100 

3.55 

0.638 

1.91 

-0.757 


TABLE XI: Contribution to the cross section pp —> uuHH + X from the subprocess uu —> uuHH. Again all cross sections 
are in fb and the errors in cos(a/) are less than or equal to 0.01. 


Ti(TeV) 

OB 

ctt 

otot 

COS (oli) 

8 

2.10xl0~ 2 

5.80X10” 3 

7.80 x 10 -3 

-0.861 

13 

6.65x 10^ 2 

1.69x 10^ 2 

2.69x 10^ 2 

-0.843 

14 

7.81 xl0~ 2 

1.96x 10“ 2 

3.20x 10~ 2 

-0.840 

33 

0.404 

8.95x 10^ 2 

0.186 

-0.809 

60 

1.08 

0.224 

0.531 

-0.786 

100 

2.36 

0.464 

1.20 

-0.776 


TABLE XII: Contribution to the cross section pp —> ddHH + X from the subprocess dd —>• ddHH. 


3.4. HH production from qq —> qqHH 


The contributions from the subprocesses uu uuHHX and dd ddHH are shown in Tables XIV and XV. 
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y/s (TeV) 

&B 

<7T 

<jtot 

cos(au) 

8 

0.404 

0.128 

0.141 

-0.860 

13 

1.18 

0.344 

0.452 

-0.841 

14 

1.36 

0.394 

0.530 

-0.836 

33 

6.24 

1.60 

2.76 

-0.804 

60 

15.6 

3.72 

7.32 

-0.788 

100 

32.2 

7.36 

15.8 

-0.772 


TABLE XIII: Contribution to the cross section pp —» udHH + X from the subprocess ud —» udHH. 


y/s (TeV) 

ob 

< 7t 

&TOT 

cos(aj) 

8 

1.46xl0 -2 

3.82xl0' 3 

9.58X10" 3 

-0.592 

13 

5.03X10' 2 

1.27x 10 -2 

2.92x 10 -2 

-0.669 

14 

6.00X10" 2 

1.50xl0' 2 

3.44xlO" 2 

-0.677 

33 

0.366 

8.36 x 10~ 2 

0.196 

-0.725 

60 

1.08 

0.231 

0.577 

-0.735 

100 

2.52 

0.512 

1.36 

-0.736 


TABLE XIV: Contribution to the cross section pp —> uuHH + X from the subprocess uu —» uuHH. 


y/s (TeV) 

&B 

CTT 

otot 

cos (or) 

8 

1.33x 10~ 2 

3.56x 10“ 3 

8.95xl0~ 3 

-0.575 

13 

4.74xl0~ 2 

1.21xl0~ 2 

2.79x 10~ 2 

-0.660 

14 

5.67x 10~ 2 

1.49xl0“ 2 

3.30x 10~ 2 

-0.664 

33 

0.359 

8.30x 10~ 2 

0.193 

-0.721 

60 

1.10 

0.238 

0.590 

-0.731 

100 

2.61 

0.538 

1.41 

-0.733 


TABLE XV: Contribution to the cross section pp —> ddHH + X from the subprocess dd —» ddHH. 


4. CONCLUSIONS 


As discussed in the Introduction, the interference, as parameterized by cos(aj), does not depend on the 
value of the Higgs coupling. For two Higgs production the interference between the graphs with the trilinear 
coupling and those without is usually constructive. In the exceptional cases of gg —> HH and qq —> qqHH the 
interference is strongly destructive. For all cases the approximate constancy of cos(a/) is a result of the near 
threshold behavior of the amplitudes and the large values of the parton luminosity in the threshold region. For 
the dominant subprocess, gg —> HH , the almost maximal destructive interference between the triangle and box 
amplitudes (cos(aj) = —0.90) obtained using unitarity cutting [19, 20] tends to make the total cross section 
small and thus difficult to measure. 

We note that pp — > ttHH + X has constructive interference between diagrams with trilinear Higgs couplings 
and those without, and the interference angle cos(aq) depends on the collider energy as demonstrated in Table 
X. For collider energy less than 14 TeV, quark-antiquark fusion is the dominant source of ttHH, while gluon 
fusion makes larger contribution for collider energy greater than 33 TeV. Thus parton distribution functions 
enhance the threshold effect of qq —> ttHH (y/s < 14 TeV) and gg —> ttHH (y/s > 33 TeV) in pp —> ttHH + X 
at different values of collider energy. 

For all the processes, except gg —> HH, unitarity arguments are not of use since the leading amplitudes occur 
at the tree level. Nevertheless, these partonic amplitudes lead to a similar maximal interference in the hadronic 
cross sections near threshold. As can be seen in Tables XVI and XVII for uu —t ZHH and uu —t ttHH, the 
maximal interference region is not quite as persistent as in the dominant gg —> HH case so the integration over 
the relevant parton luminosities results in a clustering of the interference terms rather than yielding a constant 
value. 

If we allow the standard model Higgs trilinear coupling to vary by a factor n (n\v H 3 in Eq. (1)), the results 
for gg —> HH arc illustrated in Fig. 2. In addition, Fig. 3 shows that cos (a/) remains constant for any k. 
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Vi (GeV) 

ds(/fe) 

a T (fb) 

^Tor(fb) 

cos (ai) 

350 

0.453x 10' 3 

0.126x 10 -3 

0.106X10" 2 

0.999 

400 

0.122X10” 1 

0.277xl0~ 2 

0.262x 10 _1 

0.964 

500 

0.312X10” 1 

0.571x 10 -2 

0.589X10” 1 

0.821 

1000 

0.335X10” 1 

0.416x 10 -2 

0.453 x 10 _1 

0.325 

2000 

0.181x 10' 1 

0.148x 10~ 2 

0.207X10” 1 

0.106 

4000 

0.794x 10~ 2 

0.421x 10“ 3 

0.849 xlO" 2 

0.0337 


TABLE XVI: To separate the effect of the distribution functions we show the partonic cross sections a(uu —> ZHH ) 
along with cos (a/) as a function of VS- 


V§ (GeV) 

°B(fb) 

a T (fb) 

CTTOT(fb) 

cos (or) 

600 

0.469x10'® 

0.126154x 10“ 7 

0.655x10'® 

1.000 

700 

0.340X10' 1 

O.lOlxlO' 2 

0.464X10' 1 

0.975 

800 

0.119 

0.343 xlO' 2 

0.160 

0.933 

900 

0.189 

0.548xl0' 2 

0.251 

0.890 

1000 

0.229 

0.676 xlO' 2 

0.302 

0.849 

1500 

0.212 

0.714xl0' 2 

0.273 

0.690 

2000 

0.146 

0.552xl0' 2 

0.185 

0.588 


TABLE XVII: Same as Table XVI for a(uu ->■ ttHH). 


A summary of the results for cos(aj) given in the Tables is shown in Fig.4 which illustrates again that, while 
the value of the interference is very different for different processes, there is very little variation with energy 
within a given process. 



FIG. 2: The dependence of the pp —> HH + X cross 
section from the subprocess gg —> HH as a function 
of the ratio of the trilinear coupling to the Standard 
model coupling is illustrated for yfs = 14 TeV and the 
2.27 if-factor from Eq. (3). 
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FIG. 3: The invariance of cos(au) for changes in the 
trilinear couplings is shown for the particular case of 
y/s = 14 TeV and the 2.27 AT-factor from Eq. (3). 
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FIG. 4: The values of cos(o/ ) are plotted. Note that each partonic initial state has been integrated over the appropriate parton 
distribution functions for pp collisions. For each process the different lines are for different energies ranging from 8 TeV for the line 
with the largest magnitude to 100 TeV for the line with the smallest magnitude. 
























